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Abstract 

We examine whether observable majoron emission in double beta de- 
cay can be compatible with the big-bang nucleosynthesis (NS) and the ob- 
served neutrino flux from SN1987A. It is found that the NS upper bound 
on 4 He abundance implies that the majoron-neutrino Yukawa coupling 
constant g < 9 x 10~ 6 and its maximal value is allowed only when the 
scalar quartic coupling constant A is extremely small, A < 100g 2 . It is 
also observed that, although quite less restrictive, SN1987A also provides 
independent constraints on coupling constants. 



Majoron is the massless Goldstone boson associated with spontaneous lepton 
number violation ffl. If exists, it would lead to many interesting phenomenological 
consequences. Amongst them, one particularly interesting phenomenon is "majoron 
emitting neutrinoless double beta decay" {(3(3 j) @. Recently it was pointed out that 
the potential anomaly in the double beta decay spectra of several elements may be 
explained by (3(3j |§. The desired value of the Yukawa coupling constant is roughly 
g ee ~ 10~ 4 . Regardless of this observation, if g ee is not very small, e.g. not less than 
10~ 5 , one may be able to observe (3(3 j in the near future. As was shown in ref. ||, 
it is not difficult to construct majoron models which provides such a value of g ee 
while satisfying all known experimental constraints. However in view of the strong 
cosmological and astrophysical implications of majoron 0, it is desirable to examine 
whether observable (3(3j can be compatible with the big-bang cosmology and also with 
astrophysical observations. In this paper, we consider possible constraints from the 
big-bang nucleosynthesis and the supernova SN1987A on majoron models in which 
observable (3(3j is possible. 

If the lepton number symmetry L is spontaneously broken at an energy scale v, 
majoron- neutrino Yukawa coupling matrix g a p (a,(3=e, fx, r) is related to the neutrino 
mass matrix m a p as g a p = rriap/v. Current data on z/-less (3(3 decay implies m ee < 1 
eV and g ee < few x 10~ 4 with uncertainties arising from nuclear matrix elements ||. 
For (3(3 j to be observable in the near future, we may need g ee > 10~ 5 . This implies 
that L-breaking scale is very small compared to the Fermi scale, 

v < 100 keV, (1) 

leading to a fine tuning problem in general ||. Here we will not concern this theo- 
retical difficulty. We rather concentrate on models with such a low L-breaking scale 
to see whether observable (3(3 j can be compatible with the standard nucleosynthesis 
model and the observed neutrino flux from SN1987A ||. 
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In majoron models adopting low energy L- violation (1) to provide observable (3f3j, 
low energy majoron- neutrino interactions can be described by the effective lagrangian: 

Ceff = v a iTdv a + \d^x\ 2 - 4=(^X*^>/3 + h.c.) 

-A(xx*-y) 2 + - (2) 

where \ is a (mostly) gauge singlet Higgs field carrying the lepton number two, and 
the ellipsis denotes generic nonrenormalizable interactions suppressed by the powers 
of cutoff scale A which is usually taken to be the Fermi scale. The majoron field J 
appears in \ as 

X=-j=(v + P + iJ), (3) 
and the mass eigenstates neutrino z/j = J2a Ui a v a (i = 1, 2, 3) has the Yukawa inter- 
action with J and p whose coupling constant is given by Qi — \ J2 a /3 U* a U*pg a p\ for a 
unitary mixing matrix Ui a . 

Let us now consider possible constraints on majoron-neutrino interactions from 
the big-bang nucleosynthesis (NS) 0. As is well known, in the standard model of 
NS, the energy density of exotic particles at the NS epoch (T/vs — 1 MeV) is severely 
constrained |TIJ. In the case that L is restored and also there is no sizable lepton 
number excess, the observed 4 He abundance implies 

Y(T NS ) = p x (T NS )/p u (T NS ) < 0.3, (4) 

where p x and p v denote the energy density of x an d a single species of left-handed 
neutrino respectively jlli . Since v < T^s, before the NS the ratio of the Hubble 
expansion rate H to the v-\ interaction rate behaves as T- mt /H ~ T~ l . Then x's 
would be rare at high temperature, but eventually enter into a thermal equilibrium 
with neutrinos at some temperature T eq at which r; nt ~ H. For the NS constraint 
(4) to be satisfied, we need T eq < T NS - Here we will directly evaluate p x and apply 
the constraint (4), rather than using more naive condition T eq < T/vs- 
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To evaluate p x , we first need to know whether L is restored around the NS epoch. 
The effective mass of x i n the early universe is given by [ O 



d 3 k 
(2tt) 3 E' 

~Xv 2 + 4Xn x (1/E X ) + tr(g/)T 2 /12, (5) 



"4 = -Xv 2 + J-^m x + 2tr(ggl)U) 



where fx (X = \i z/ ) denotes the phase space distribution function of X whose 
number density is defined as rix = J d 3 k fx/(2ii) 3 , and neutrinos are assumed to be 
in thermal equilibrium. To proceed, let us set 

n x (1/E X ) = £Y"n u (1/E U ) = ^T 2 /24, (6) 

where Y = p x j ' p v . Clearly the average energy of x do not exceed that of v and thus 
n x /n v > Y, implying that > Y. Then for v < 100 keV and Y saturating the 
NS bound (4), which is the most interesting case for us, L is restored around the NS 
epoch regardless of the value of tr(ggt) fT3|. 



It may be necessary to further discuss on the parameters £ and io. Self interac- 
tions among x's do not change p x , but can increase n x through the processes like 
XX XXXX* ■ The values of £ and us depend on the strength of such "number 
changing self interaction process" (NCP). In our case, NCP's can occur through the 
quartic coupling Xx 2 X* 2 (loop effects) or through nonrenormalizable interactions like 
K X 3 X* 3 7 ' X 2 '. If the NCP rate is weaker than the expansion rate, the average energy of 
X is roughly equal to that of z/, and then £ ~ u ~ 1. In the opposite case, x' s would 
achieve a thermal distribution, giving £ — 2 and uo ~ 1/2. In the subsequent analysis, 
we will simply set £ ~ uj ~ 1 since as we have argued Y < ZY^ and then this choice 
gives more conservative result for the x-production. 

Clearly x's are mainly produced by the heaviest mass eigenstate neutrino v which 
has the largest Yukawa coupling g = max (gi). For interaction terms in (2), the 



processes which can dominantly produce x are as follows: (A) vv — > xj (B) vv — > XX*> 
(C) 1/1/ — > xxx* an d vv — > XX*X*- Before x' s enter into an equilibrium with z/'s, 
particularly when the NS constraint (4) is satisfied, we can safely ignore the inverse 
processes annihilating x's. Then the Boltzmann equation describing the evolution of 
p x before the onset of the NS is given by @ 

^p x + AHp x ^J2l[dX](2 7 r)^(p i -p f )E I \M I \'f J (7) 

where [dX] = Y\ a d 3 k a /(2ir 3 )2E a (a runs over particle species participating in the 
process), Ej (I = A, B, C) is the total energy of x's in the final state of the process J, 
\Mj\ 2 is the amplitude squared including the symmetry factor for identical particles, 
and finally f = Tiifi f° r the phase space distribution functions of particles in the 
initial state. After a straightforward computation, this Boltzmann equation can be 
cast into the following form: 

- H^Y = 6.7 x 10-y A[(Z) A + (1.6^ \n(5/XZ)) B + (l(T 2 A)c], (8) 

where Z = Y + (g 2 /2X) — 6(v/T) 2 . Here the subscripts in the brackets denote the 
contributing process. Note that the rate of (A) is proportional to m 2 s oc Z, while 
that of (B) includes the factor ln(mg ff ). Of course the above equation is valid only 
when L is restored, viz Z > 0. At any rate, it indicates that if A > g 2 lng~ 2 , x's are 
mainly produced around the NS epoch by the inverse decay process (A) whose rate 
is dominated by the background matter induced piece which is proportional to Y. 
However for A < g 2 lng~ 2 , the process (B) dominates over other processes. If A ^> g, 
the process (C) can be important at the very early stage of Y 1. 

The Boltzmann equation (8) can be integrated to obtain F(Tjvs). Applying the 
NS constraint (4), we then find 

A# 2 < 7.2 x 1(T 19 ln(l + e\/g 2 ), (9) 
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where e = 30/(A 2 /# 2 - 1601n(^ 2 /5 + 0.1A)). This gives 

g<9x 10~ 6 R (10) 

where R ~ 1 for A < 100g 2 , while R ~ (r _1 lnr) 1/4 for r = A/lOOg 2 > 1. Note that 
in the case of A > 100g 2 , the ^-production is dominated by the process (A) whose 
rate is enhanced by the factor X/(g 2 In g~ 2 ) compared to the rate of the process (B). 
This is the reason why we get a stronger bound on g in this case. 

The implication of our NS result for (313 j is clear. Since g is the majoron Yukawa 
coupling constant of the heaviest mass eigenstate neutrino, g ee < g and thus any 
upper limit on g applies also to g ee . Then taking into account possible uncertainties 
of our analysis, we can conclude that g ee ~ 10~ 5 for which observable j3j3j is barely 
possible can be consistent with the standard NS model, but only under very unlikely 
conditions that the quartic coupling constant A is extremely small (< 100g 2 ) and 
also the majoron Yukawa couplings with and v T do not exceed that with v e . This 
conclusion is valid as long as v is small enough, e.g. less than few hundreds keV, for 
L to be restored around the NS epoch. 

Although widely accepted and very natural, the standard NS model is not a 
unique model explaining observed cosmological data. There may be other successful 
models in which the constraint (4) is not valid any more |7]], implying that the NS 
limit on g does not have a strictly firm foundation. In this regard, it is still worthwhile 
to consider other implications of observable f3(3j with g ee > 10~ 5 , e.g. for supernova 
dynamics. It has already been studied how majoron-like particles can affect the 



explosion and the subsequent cooling of supernovae [14, pj|, but the results of these 



studies do not show any meaningful implication for (3(3 'j. Just after the observation of 



SN1987A, Kolb and Turner [16] derived a limit on the interactions between supernova 



neutrinos and "cosmic background majorons" (CBM). Since the observed neutrino 
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pulse from SN1987A is that of ui, the mass eigenstate neutrino which is mostly u e , the 
relevant Yukawa coupling constant here is g\. Although in principle the parameter 
g ee describing (3(3 j can be significantly different from g\ (note that g ee = J2i Uf e gi), it 
is somewhat natural to assume that neutrino mixing is small enough for gi to be close 
to g ee . In this regard, the consideration of CBM's may provide a useful information 
on (3(3j. 

If V\S from the supernova are scattered off by CBM's, there would be a substantial 
decrease in the average neutrino energy, leading to an effective loss of detectable flux. 
The observed data || indicates that the mean free path / of V\ through the CBM's 
should be comparable to or greater than the distance to the supernova, viz 

I > 1.7 x 10 23 cm. (11) 



Using this, it was found in ref. [16 that g% < 10 3 . However as we will see, one can 
in fact obtain a significantly stronger limit unless A is significantly less than 10~ 3 . 

Clearly for g ee in the range of observable (3(3 j, i.e. g ee > 10~ 5 , there was a period 
in the early universe when p and J were at thermal equilibrium with neutrinos. Later 
that, the relic majoron number density can be increased by the decays of p and v 
and also by the v annihilations. With these observations, we will set the majoron 
temperature at present as Tj ~ 1.9 K, which is a somewhat conservative choice. Then 
the inverse mean free path of V\ propagating through CBM's is given by 



r 1 



V2 
4vr 2 



poo 1 fl 

1* dww 2 ^—- dz{l-zf' 2 o{s), (12) 

J 6 \. J — 1 



where a(s) is the total cross section for the reaction v\ J — > v\ J with the total energy- 
momentum squared s = 2ETjw(l — z) + m\. Here E (~ 10 MeV) and mi denote 
the energy and the mass of the incident v\. There are three diagrams responsible for 
the reaction V\J — ► v\ J. Two of them are the Compton-type ones, while the rest one 
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involves the p-exchange. The resulting cross section can be written as 




[F(s) + y- 1 G(s)} 



(13) 



where 



G(s) 



5 2 2(1 + 2x) 2 ln(l + x) A+(l-x)\ny 
2 x x 2 {l + x) 1 + x 

ln(l + x) 1 



x 2 x(l + x) ' 



for y = m\/s and x = s/m 2 = s/2Xv 



2 . Here F(s) represents the contribution from 



the Compton-type diagrams, y 1 G(s) is from the p-exchange diagram. 

If A is so small, e.g. A <C gi^/v, that the p-exchange diagram can be ignored, 
we have a ~ gf ln(s/m 2 )/lQirs. Then applying the supernova constraint (11) yields 
gi < 10~ 3 as was obtained in ref ||16|| . However in the other case of A 3> gi-^s/v, the 
cross section is largely enhanced by the p-exchange diagram, allowing us to obtain 
a more stringent limit on g^. Our results summarized in fig. 1 shows that, unless 
A < 10 -3 , some portion of the parameter region giving gi < 3 x 10 -4 is ruled out by 
the supernova constraint (11). As was mentioned, if neutrino mixing is small, so that 
g ee ~ gi, this supernova result can be directly applied to (5f3j. 

To conclude, we find that the standard nucleosynthesis model strongly constrains 
the coupling constants in majoron models in which observable "majoron emitting 
z/-less double beta decay" (f3(3j) is possible through a spontaneous lepton number 
violation below few hundreds keV. Our results directly applies to the Yukawa coupling 
constant g ee governing /3/3j, yielding g ee < 9 x 10 -6 . Here the maximal value 9 x 10 -6 is 
allowed only when (i) the scalar quartic coupling is extremely weak, roughly A < 10 -8 , 
and (ii) g ee is rather close to the majoron Yukawa coupling constant of the heaviest 
neutrino species. We also find that the consideration of the scatterings between cosmic 
background majorons and supernova neutrinos provides a constraint on the majoron 
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Yukawa coupling constant g\ of the mass eigenstate neutrino which is mostly v e . 
Compared to the nucleosynthesis constraint, this supernova constraint is much less 
restrictive, but is independent of the validity of the standard nucleosynthesis model. 
If gi ~ g ee and A is not less than 10~ 3 , which is somewhat plausible, some part of the 
parameter region giving observable [3f3j is excluded by the supernova data. 

This work is supported in part by KOSEF through CTP at Seoul National University. 
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Figure 1: Supernova constraint depicted on the plane of the L breaking scale v and 
the electron neutrino mass mi. The solid curves from bottom to top are obtained 
by applying the condition (11) for A = 1CT 3 , 1CT 2 , 1CT 1 , and 1 respectively, and the 
regions below them are excluded. The two dotted lines correspond to g x = 1CT 3 and 
1CT 5 . 
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